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A variant of compression optical coherence elastography for mapping relative tissue stiffness is reported. Unlike 
conventionally discussed displacement-based (DB) elastorgaphy, in which the decrease in the cross-correlation 
is a negative factor causing errors in mapping displacement and strain fields, we propose to intentionally 
use the difference in the correlation stability (CS) for deformed tissue regions with different stiffness. We 
compare the parameter ranges (in terms of noise-to-signal ratio and strain) in which the conventional DB- 
and CS-approaches are operable. It is shown that the CS approach has such advantages as significantly 
wider oper ability region in terms of strain and is more tolerant to noises. This is favorable for freehand 
implementation of this approach. Examples of simulated and real CS-based elastographic OCT images are 
given. © 2012 Optical Society of America 

OCIS codes: (110.4500) Optical coherence tomography; (170.6935) Tissue characterization. 



The problem of elasticity imaging in optical coherence 
tomography (OCT) has been attracting much attention 
since the end of 90ies [1-9]. However, the elastography 
regime is not yet realized in commercial OCT scan- 
ners, whereas in medical ultrasonics combined conven- 
tional scans and shear-elasticity imaging are already im- 
plemented in several platforms (e.g., Siemens, Hitachi, 
Ultrasonix [10]). Although some ultrasonic systems use 
excitation of shear waves and measuring their velocities 
(e.g., Fibroscan and Supersonic scanners) [10], mostly 
the elasticity imaging is based on special processing of 
conventional ultrasonic images. Since the first works on 
OCT-elastography [1], similar processing of OCT images 
is often discussed to reconstruct displacements produced 
by additional compression/shearing of the tissue. Sub- 
sequent differentiation of the displacement field can be 
used to determine the elastic-strain field and thus to es- 
timate the shear- modulus distribution [1,2]. 

Despite the apparent simplicity of the idea, the imple- 
mentation of this approach faced significant difficulties 
and required highly controllable laboratory conditions 
in experimental demonstrations (e.g. [1,2]). Indeed, nu- 
merical differentiation is rather error-sensitive and there- 
fore requires accurate reconstruction of the displacement 
field. To fulfill this requirement, the deformations should 
not be too small. Otherwise, the pixelized structure of 
OCT images does not ensure acceptable accuracy even 
with application of various super-resolution smoothing 
techniques [10,11]. On the other hand, the tissue strains 
should not be too large, because distortions of the scat- 
terer patterns in the deformed OCT images reduce the 
accuracy of the cross-correlation (i.e., produce "decorre- 
lation noise" [11]). Inevitable other noises in real images 
introduce additional complications. In what follows we 
show that the above-mentioned factors limit the accept- 
able upper strains to several percentages. This limitation 



significantly complicates realization of the most practi- 
cally interesting freehand mode of elastography in OCT. 

Here, we report an alternative variant of mapping the 
relative stiffness using comparison of OCT images ob- 
tained under different degrees of straining. We call this 
variant the correlation-stability (CS) approach, bearing 
in mind that stiffer regions experience smaller distortions 
and demonstrate higher cross-correlation (which was also 
mentioned for ultrasonic scans [12] ). After presenting 
numerical simulations and comparison with the conven- 
tional DB approach we give some in vivo examples of CS 
maps of the relative elasticity distributions obtained in 
freehand mode using a spectral-domain OCT scanner. 

For the CS mapping, the reduction of cross-correlation 
between the images of the deformed tissue is an infor- 
mative factor, rather than "decorrelation noise" reduc- 
ing the mapping accuracy. There is a similarity with the 
nonlinear-acoustic approach to detection of cracks (see, 
e.g., [13]), in which nonlinear distortions of the sound- 
ing field produced by increased nonlinearity of the de- 
fects are intentionally used as a signature of their pres- 
ence. Basically the procedure of cross-correlation be- 
tween images using a moving (usually rectangular) win- 
dow mi x vfi2 in size is similar to that used in refs. [1,2]: 
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where i = l..mi, j = l..m2, the window is moved by n 
pixels axially and k laterally, jap^R are the mean values in 
mi x rri2 areas on images P and R. For ideally coinciding 
areas (7(0,0) = 1 and tends to zero for uncorrelated 
ones. For real partially distorted and displaced areas, the 
position (n*,fc*) is found where the correlation reaches 
maximum (7* and the field C*(i,j) is plotted as CS-map. 
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To give an example of the CS approach let us consider 
a simulated image of 400x200 pixels in size (typical of 
OCT images). The "scatterers" are formed by putting 
initially a random value in each pixel. Then Fourier fil- 
tering is applied to smooth the image and obtain inhomo- 
geneities with correlation properties similar to that typ- 
ical of real OCT images. Namely, the filtering has to en- 
sure the same as in real OCT images level of correlation 
between non-overlapped correlation windows containing 
independent scatterers. This background correlation C bg 
is determined by the number q\ x q2 of independent 
scatterers within the correlation window (qi,2 < ^1,2)- 
It can be shown [14] that C bg = (tt /2) 1 / 2 / '(qiq 2 ) i/2 . 
In real OCT images, e.g., of a human skin, typically 
C bg ~ 0.3 - 0.2 for window sizes mi 5 2 ~ 20 — 40 (similar 
windows were used in [2]). Such values of C bg correspond 
to sizes of independent scatterers about several pixels, 
which was ensured for the simulated images. 

The next step is simulation of strained- tissue images. 
In real experiments the surface of the OCT sensor acts 
as a rigid piston producing strain fields combining axial 
and lateral components. Here, for illustration we con- 
sider purely shear lateral deformation. This allows one 
to use a simple analytical solution even for inhomoge- 
neous tissue containing a layer with a contrasting shear 
modulus /i2 = b • /i, where \i = \i\ = \±<$ is the modulus 
of the surrounding upper and lower layers, and b is the 
contrast. If the total sample thickness is I = l\ + I2 + Z3, 
where I2 is the contrasting-layer thickness (see Fig. la) 
and the lateral displacement of the tissue surface is uo, 
then the lateral displacements u(z) inside the sample are 
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and the mean shear strain is 
e = uq/1. The patterns of the scatterers in the sheared 
sample can be recalculated using Eqs. (2). Continuous 
subpixel displacements can also be correctly found by 
using Eqs. (2) and applying forward and inverse Fourier 
transforms of the pattern in combination with the theo- 
rem about the Fourier spectrum of shifted functions. 

Figure la shows the initial and deformed patterns of 
simulated scatterers in which the contrasting-stiffness 
layer is not yet visible. An example of CS mapping of this 
layer is shown in Fig. lb, where the stiffer layer is clearly 
seen due to increased correlation compared with the sur- 
rounding softer layers that experience stronger distor- 
tions. Finally, Fig. lc shows the reconstructed displace- 
ments in a vertical slice using the conventional cross- 
correlation approach (with parabolic smoothing of cor- 
relation peaks to obtain sub-pixel resolution). In Fig. lc 
the stiffer layer looks as a butterfly-like deviation of the 
u(z) function (solid line) from the dashed line u m (z) 
corresponding to the mean shearing of the entire sam- 
ple. Conventionally [1], the so- found displacement field 
is supposed to be differentiated to obtain the strain field 
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Fig. 1. Simulated example of mapping the stiffer layer 
with contrast 6 = 3 and thickness I2 = 20 pixels. Plot 
(a) shows the pattern of scatterers with C bg ~ 0.3 for 
e = 0% and e = 50%, (b) is the CS image of the stiffer 
layer obtained by cross-correlating the unstrained pat- 
tern with deformed one (see plot (a)). Plot (c) shows the 
mean displacement u m (z) (dashed line), reconstructed 
lateral displacement u{z) (solid line) along one of verti- 
cal paths, and its derivative averaged over all such paths 
(dotted line) found by the conventional cross-correlation 
technique for the mean shear strain 6%. The correlation- 
window size is 20x20 pixels. The images contain no ad- 
ditional noise and the errors are exclusively due to de- 
formation of the scatterer pattern. 



(dotted line shown in Fig. lc in arbitrary units), in which 
regions with smaller strains correspond to larger stiff- 
ness. Such a differentiation is known to introduce ad- 
ditional errors (even if there are no other noises than 
the deformation-produced "decorrelation noise" [11]). In 
Fig. lc the differentiated dotted line even after averaging 
over all vertical cuts demonstrates pronounced irregular- 
ities which are exclusively due to strain-produced distor- 
tions of the scatterer pattern. Hereafter, all distortions of 
the image unrelated to the stress-produced deformation 
of the scatterer pattern we call noise. 

Now we compare the operability regions for the pro- 
posed CS mapping and conventional DB elasticity map- 
ping in terms of strain and tolerance to additional noise. 
In our comparison we even adopt the form more favor- 
able for the conventional approach, i.e., avoiding error- 
sensitive differentiation and bearing in mind that the dis- 
placement field itself can already allow one to detect the 
presence of stiffer regions (which is also mentioned in [1]). 
Namely, we use the presence of the "butterfly" shown 
in Fig. lc as a signature of the presence of the stiffer 
layer. The mean slope (see the dashed line in Fig. lc) 
can be subtracted for clarity to obtain a horizontally- 
oriented "butterfly" as shown in Fig. 2a. The informa- 
tive "butterfly" with amplitude A\> can be masked by 
various irregularities characterized by the variance a u 
due to the presence of noises and decorrelation distor- 
tions of the scatterer pattern in OCT images. As the 
threshold condition we adopt the often used in statistics 
"three-sigma criterion", i.e., A^ > 3a u . The variance o u 
of the irregularities is estimated for each given strain 
by gradually increasing the amplitude of the generated 
additional noise and then comparing a u with A^. For 
the CS mapping, we make a similar comparison between 
the cross-correlation level C sU f f for the stiffer layer and 
(j so ft £ or ^Yiq surrounding softer areas. These two lev- 
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Fig. 2. Elucidation to the determining the operability 
areas for the conventional DB- and CS- mapping. Plot 
(a) shows a bunch of vertical slices for u{z) — u m (z) with 
the mean "butterfly" of amplitude corresponding to 
the stiffer layer and irregularities in the reconstructed 
displacements with variance a u ; (b) CS map with regions 
having difference in the mean correlation < C sU ^ — 
(jsoft > anc j backgound variance cr h c 9 orr ] (c) operability 
regions of the two methods in terms of strain and noise- 
to-signal ratio N/S defined as the ratio of square roots 
of the noise and signal energies (found after subtraction 
of the respective mean values ). 

els of cross-correlation are characterized by their respec- 
tive average values and the background cross-correlation 
variance (J b c 9 orr (which is due to irregularities in the cross- 
correlation coefficient that are always present even in 
the absence of the stiffer layer). For a given strain, we 
vary the noise level and again determine the threshold 
difference < C stiff - C soft >= 3a b c g orr . The so-found 
curves delimiting the operability regions for the DB and 
CS mapping are shown in Fig. 2c. In agreement with 
other authors [3,11] we found that the conventional DB 
mapping works for fairly small (below ~ 10%) strains, 
whereas the CS mapping works in a much wider range of 
larger strains (~ 20 — 150%), which is favorable for free- 
hand operation. Besides, CS mapping is noticeably more 
tolerant to noises. For the demonstrations, we used a re- 





Fig. 3. Two in vivo examples of CS maps where red color 
shows the stiffer regions. Panel (a) shows a conventional 
OCT image where a hair root is hardly seen and (b) is 
the corresponding CS map for correlation- window size 
20 x 20 pixels. Panels (c) and (d) show OCT images of 
another bulb superposed with the elastographic CS maps 
with different resolution obtained using windows 40 x 40 
and 20 x 20 pixels, respectively. 



cently developed OCT spectral-domain scanner based on 
techniques described in [15,16]. It acquires 480x289 pixel 
images with a rate of 21fps, axial resolution 10 /im and 
lateral one 20 /im. The fiber-optic OCT probe is 2 mm 
in diameter and acts as a piston producing deformation 
of the tissue during freehand operation. The images are 
fed in a PC in real time and after filtering, amplitude 
normalization, and averaging of a few frames are cross- 
correlated with a reference one. For in vivo examples we 
have chosen human cheek skin where hair roots serve 
as stiffer inclusions. The resulting CS maps are shown 
in Fig. 3, where the stiffer regions are clearly seen. For 
these examples, the average strain created in the tissue 
by pressing the probe is about 20 — 50% in agreement 
with the above-found operability region for CS mapping. 

Thus the performed analysis and tests demonstrate 
that the rather straightforward CS approach can readily 
be implemented in sufficiently high-speed OCT scanners 
in freehand mode to obtain relative-stiffness maps simi- 
lar to those used in medical ultrasonics [10]. 

Support of the grant 11.G34.31.0066 of the Russian 
Federation Government is acknowledged. 
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